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ABSTRACT 


This report presents a theoretical discussion of the effect of non¬ 
equilibrium quantities of free radicals on the ignition delay in the 
hydrogen-air system. A correlation was found between the ignition 
delay computed from, complete reaction kinetics and the initial rate of 
the reaction 0, + H = OH + 0 . A rate equation is derived in which a 
modified expression for the total free radical content is substituted for 
the H concentration. The implications of non-equilibrium amounts of 
free radicals on the testing and operation of high velocity combustion 
systems are discussed. 
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1.0 INTRODUCTION 


Currently, there is considerable interest in the ignition and confi- 
bustion processes which occur in supersonic flow (Refs. 1-6). In most 
cases of interest, the ignition temperatures are relatively low, typi¬ 
cally under 2000'’K as opposed to greater than 2000°K for shock tube 
experiments. Two reasons for the interest in low temperature ignition 
phenomena are (1) very high initial temperature means a high combus¬ 
tion temperature wdth attendant dissociation losses and (2) even 
moderately high static temperatures are very difficult and expensive to 
obtain in steady flow facilities. With these relatively low temperatures, 
ignition delay can occupy a sizable fraction of the length of the combus¬ 
tion chamber, and therefore, the factors which affect ignition delay be¬ 
come important. 

In a classical ignition reaction in which the reactants are relatively 
slowly heated, ignition will occur when the equilibrium free-radical 
concentration reaches a point where the rate of formation of free radi¬ 
cals by chain branching reactions exceeds their loss by recombination 
on the vessel walls. Except in the case of adiabatic compression, heat 
transfer rates prevent heating a gas fast enough so that the free-radical 
content is not in equilibrium with the parent species. In supersonic flow 
where rapid adiabatic compressions and expansions are the rule, it is 
almost impossible to perform these processes slowly enough so that 
equilibrium will be reached, and therefore, the reacting mixtures may 
depart radically from their equilibrium free-radical contents. 

The purpose of this investigation is to present the effect of non¬ 
equilibrium free-radical content on ignition properties and to determine 
the difficulties or benefits which may result from these effects. This 
investigation was made in the Research Branch of the Rocket Test Facility 
(RTF), Arnold Engineering Development Center (AEDC), Air Force 
Systems Command (AFSC). 


2.0 DISCUSSION 


In the hydrogen-air system, the length of the ignition delay is con¬ 
trolled by the rate of increase in free radicals through the four following 
reactions (Ref. 7); 


0,' + 

H = OH +0 

(1) 

20H = H,0 + 0 

(3) 

0 + 

H, = OH + H 

(2) 

OH + = H,0 + H 

(4) 


Manuscript received October 1964. 
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No I’eactions with niii-o^en are considerecl in this work. The rate of 
formation of a particular radical (.A / from reactants (H) and (C) through 
a particular reaction (i) is 



I’Ai 


_ - K i / R T 

lir 


(X FI) (X c) 


The rate of the reaction will depend on the temperature and reactant 
concentrations; its sensitivity to temperature w'ill depend on the activa¬ 
tion energy (Ei). 


- In a supersonic test facility, it is possible to have the free-radical 
content of the gas stream much larger than the value at equilibrium at 
the local static temperature. This would occur when the gas had been 
heated and rapidly expanded through a nozzle. The free-radical con¬ 
centrations could be nearly in equilibrium with conditions at the nozzle 
throat. It could also occur to an even greater degree in systems heated 
by an electrical discharge where greater than equilibrium amounts of 
dissociation products could occur in the arc and remain frozen during 
the expansion process. 

If the concentration of free radicals is high enough to cause the 
initial rate of reaction to exceed some critical value, the reaction will 
proceed exponentially to a local equilih'rium condition. Reactions (!)• and 
(2) feed each other with free radicals and (3) and (4) produce m.ore reactants 
for (1) and (2) from the OFI produced b}^ (1) and (2). The set of reactions will 
proceed until an equilibrium is reached between the reactions as shown and 
the reverse reactions or until the free-radical concentrations become so 
high that the recombination reactions become important. The critical 
value of the reaction rate depends on several factors which depend on the 
recombination rates of the radicals involved and on the time available for 
ignition. Th-^ time factor is important in flow S3^stems since, at 10,000 ft/ 
sec with an ignition delay of 0. 01 sec, ignition would occur 100 ft from the 
point of initiation and for practical purposes might as v/ell not occur at all. 

The rates of the four principal reactions in the forward direction are: 


= 4 


10 ' 


RT 


18000 I 

(Xo,)(Xh) ' 


6000 


= 3 X 

10“ -L 

RT 

KT 

e 

(Xiij) (Xq) 



2500 


Tj = 3 X 


e 

(XoH )* 


2500 


r, 3 X 

if 

'‘UT 

e 

(Xoil) (Xiij) 


Ref. 8 
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In the low temperature ranj^e (< 2000°K) reaction (1) will be, by far, 
the slowest because of the large activation energy and will control the 
overall rate of the reactions since it provides the only mechanism by 
which Oj may react. For example, at 1000°K and equivalent reactant 
concentrations and pressure, reaction (2) will occur 300 times as fast 
as reaction (1). 

The rate of formation of 0 by reaction (1) is proportional to Xji . 
However, when the temperature is low, it may be assumed that in the 
very early parts of the reaction, reactions (2) and (4) will proceed nearly 
to completion before (1) will start. Reaction (3) is neglected since it is 
normally several orders of magnitude slower than (4) by the ratio of Xqh 
to Xhj since k, = k., . 

Since, in the very early part of the ignition delay period, XHj 
changes only slightly and may be assumed to be constant, an analytical 
solution for the rates of reactions (1) and (2) may be obtained from the 
following equations: 

[k, (XhJ(Xo) + K4CXoh)(Xh,)] 

[k. (Xh,)(Xo) - K,(Xoii)(Xh,)] 


and 


(Xoh)+(Xh) = 2(Xo), - 2(Xo) + (Xoh)„ + (Xh)o 


which is: 


K; P (XH, )r 

(Xo) = (Xo)„ e 


( XOH ) 


K.(Xo)^ 
K, - K, 


Kj P (XHj)r 
RT 


K,P (Xh,)' 
RT 


'(XoH). 


K^P (XHj)i 
RT 


(Xh) = 2(Xo)„ - 2(Xo) + (Xoi!)o - (Xoh) + (Xh)„ 


3 



AEDC-TDR-64.241 


The results of these equations for a stoichiometric mixture of II, 
and Oj -with 0 or OH at SOO^R are shown in Figs, la and b, respectively. 
They show that under low temperature conditions, the 0 and Oil will 
react to produce II in the very early part of :he ignition delay period 



Fig. 1 Free-Radical Concentrotions in the Very Early 
Part of the Ignition Delay Period 


wdth each 0 resulting in two H atoms and each OH in one H atom. If the 
reverse reactions had not been neglected, the OH would tend toward the 
equilibrium value rather than zero, and the amount of H produced would 
be less. However, it may be shown by the equilibrium constants for the 
reactions that the reactions will essentially go to completion in the for- 
w'ard direction. 

For the reaction OH + HjO—^H (reaction 4) , the equilibrium constant is 

„ (Xh,o)(Xh) 17000 at 1000°K 
^ (XohHXh,) 2 X lo‘ at 500°KJ^^®‘ 

And for the sum effect of reactions (2) and (4) 2Hi + 0 HjO + 2H , the 
equilibrium constant is 

V _ <Xh,o)(Xh)V _ 8500 at 1000°X \ 

(XhJ' (Xq) 2.2 X 10‘ at 500 °kJ 
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(Since the equilibrium constants are large, the reactions will go 
essentially to completion in the direction shown. Therefore, under con¬ 
ditions where reaction U) is slow compared to the others, the 0 and OH 
will react to form H by reactions (2), (3), and (4) , and (XiOoWill be approxi¬ 
mately equal to (Xh)^ + 2(Xo)(, + (Xoh)^ very early in the ignition delay 
period. ) 


Under the.=:e conditions, the initial rate of reaction (1) will be: 


, P Ki(Xoj)„ (XfB ) ,, “r°t' 

(rj)^ = --- where Kj = 4 x 10 e 


(Ref. 8) 


(Xfr)„ = (Xh)„ + (Xoh)„ + 2(Xoh 


The correlation between the ignition delay time (chosen as the time 
where (XH,) = 0.5 {Xh 2 )o), obtained from solutions to the General Applied 
Science Laboratory Computer program for the kinetics of the HjO, system 
(Ref. 7) and the initial rate of reaction (1), (ri)^, is shown in Fig. 2. If 
the initial conditions are known, the delay time can be predicted from the 
plot within a factor of 2 . 


• Simulated Combustion Heated High and 
o Arc Heated High Xq 
^ Equilibrium Air 
T ■ 800 to 1700°K 
P -0.35 to 1.85 atm 



Ignition Delay Times Pressure. Pt, atm-sec 


Fig. 2 Correlation between the Porameter 
(rdp and Ignition Delay Time 


If the wide range of starting conditions are considered, the corre¬ 
lation seems good. Mixtures with an excess of free radicals seem to 
have a somewhat larger delay time for the same (rjl^^ than equilibrium 
mixtures. Pressure seems to have large effect on delay, whereas the 
concentration terms seem to have a smaller effect. 
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If an ignition delay of 50 ^ sec is selected as a maximum which can 
be tolerated because of equipment size, a minimum (ril^, of 0. 2 must be 
maii]tained at 1 atm. By using this as a critf-ria, a relationship may be 
obtained for a minimum ignition temperature for any free-radical con¬ 
centration. This relationship is shown in Fig. 3 for a case where par¬ 
tially dissociated hydrogen is mixed stoichiometrically with air. 



400 W) 800 1000 1200 

T. °K 


Fig. 3 Effect of Non-Equilibrium H 2 Dissociation 
on Ignition Temperature 

Another example of the effect of non-equilibrium free-radical con¬ 
centrations on ignition is shown in Fig. 4. In this case, air was assumed 



Xu, mole/mole 


Fig. 4 Effect of Arc Heating Air on the Ignition 
of H 2 - Air Mixtures 
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to liav'e lieen arc heated and expanded tlirough an i.sentropic nozzle with 
a frozen coni[)osition equal to the ecjuilil^rium composition at 100 atm 
and '•!ri00°K. It was then mixed with llj at 280‘’I\ in such a manner that 
the static entlialpy of the mixture was equal to the static enthalpy of each 
component times its respective mole fraction (no velocity change on 
mixing). With the non-equilibrium a.r, ignition would occur in all mix¬ 
tures up to about 75 percent llj, whereas with equilibrium air at the 
same static temperature, there would probably not be ignition above 
about 20 percent Hj. These results would be modified somewhat in a 
real mixing zone where the velocity and, therefore, the static enthalpy 
would change with the mixture ratio. 

There are two factors which cause a large uncertainty in the con¬ 
clusion obtained from the above results. First, there is no assurance 
that the ignition criteria of (ri), is valid in as extreme a case as SOO^K 
and hydrogen dissociation of 2 percent as .showm in Fig. 3. Second, the 
rate constants and reaction mechanisms used are accepted as valid at 
normal combustion temperatures (Ref. 8), but there is no assurance 
that they apply at 1000°K and below. However, it is reasonable to assume 
that a large effect of non-equilibrium free radicals on ignition exists and 
that ignition can occur much below the normally accepted temperatures 
when the free-radical content is high. 

3.0 CONCLUDING REMARKS 


It has been shown that the presence of non-equilibrium quantities of 
free radicals in a mixture of Hj and air can greatly affect the rate of the 
reactions involving these species and, therefore, the time required for 
ignition. There are a number of implications which may be drawm from 
this conclusion with regard to the design and testing of propulsion systems 
operating with high velocity in the combustor. 

1. In a ground test facility where the air must be heated to a high 
total temperature, the gas composition will freeze during the 
expansion at some temperature above the test static temperature, 
and the excess of free radicals will speed up the ignition delay 
reactions which occurs on mixing the air with fuel. 

2. Conversely, in a free-flight system, the concentration of free 
radicals will lag the equilibrium value during the compression 
in the diffuser. This may result in a substantially reduced rate 
of chemical reaction when the mixing with Hj occurs. 
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3. In many research investigations where the objective is to in¬ 
vestigate a process such as the effects of heat release on 
aerodynamics or mixing rates, it may not be important if the 
ignition chemistry is somewhat modified. Non-equilibrium 
amounts of free radicals could be provided by partial vitiation 
or seeding with arc heated gas and used to circumvent the very 
high total temperatures necessary to get sufficient static tem¬ 
perature to cause reaction in high velocity flow. 

4. It might be feasible to use a device such as a high frequency, 
high voltage arc to partly dissociate either or both reactants 
without much heating so that experiments could be performed 
at much lower temperatures than v/ould be possible with ther¬ 
mally heated systems. 

5. Finally, one must very carefully specify departures from equi¬ 
librium when reporting temperatures at which ignition occurs in 
a supersonic flow system. 
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